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We suggest to carry out lattice calculations of current correlators in position space, sandwiched 
between the vacuum and a hadron state (e.g. pion), in order to access hadronic light-cone distribution 
amplitudes (DAs). In this way the renormalization problem for composite lattice operators is avoided 
altogether, and the connection to the DA is done using perturbation theory in the continuum. As 
an example, the correlation function of two electromagnetic currents is calculated to the next-to- 
next-to-leading order accuracy in perturbation theory and including the twist-4 corrections. We 
argue that this strategy is fully competitive with direct lattice measurements of the moments of the 
DA, defined as matrix elements of local operators, and offers new insight in the space-time picture 
of hard exclusive reactions. 



PACS numbers: 12.38.-t, 14.20.Dh; 13.40.Gp 

I. INTRODUCTION 

Hadron light-cone distribution amplitudes (DAs) 
present the principal nonperturbative input in the 
pQCD description of hard exclusive reactions and 
are to a large extent complementary to the usual par- 
ton distributions. The existing information on DAs 
is, however, very limited. The main reason for this 
is that the DAs arc much more difficult to access ex- 
perimentally: for realistic momentum transfers, the 
contributions of interest are often contaminated by 
large nonperturbative corrections coming from the 
end-point regions in the quark momentum fraction. 

The pion leading twist DA is the simplest one and 
has attracted most attention. There is increasing ev- 
idence that at the scale of the order of 1 GeV this 
DA differs considerably from its asymptotic form. 
In particular, QCD sum rule estimates [1, 2], lat- 
tice calculations [3] and the analysis of the experi- 
mental data on the transition form factor 7* 777 
[1, 4, 5] are consistent with the positive value of the 
second Gegenbauer moment of the pion DA which 
is roughly a factor two below the original estimate 
by Chcrnyak and Zhitnitsky [6]. Beyond the second 
moment, very little is known. The analysis of the 
7* — !■ TTj form factor in Ref. [1, 4] indicates a neg- 
ative value for the fourth Gegenbauer moment, but 
the status of this result is not clear as the analysis 



has some model dependence. The lattice calcula- 
tions of the fourth and higher Gegenbauer moments 
would be very difficult because they contribute with 
a small coefficient, and because the lattice renor- 
malization of local operators with many derivatives 
becomes too cumbersome. The aim of this letter 
is to suggest an alternative approach, based on the 
lattice calculation of exclusive amplitudes in coor- 
dinate space. We will argue that the interpretation 
of such calculations in the framework of QCD fac- 
torization is equally straightforward and offers new 
insights compared to the standard momentum space 
formulation. From the lattice side, the advantage is 
that the renormalization problem for composite op- 
erators is avoided altogether, but, instead, in order 
to be sensitive to the detailed structure of the pion 
DA, one needs pion sources with large momentum, 
at least of order 2-3 GeV. 

The idea to emphasis the coordinate rather 
than momentum-space dependence of the correla- 
tion functions is by itself not new, see e.g. [7, 8]. 
Our proposal goes in the same direction as the work 
[9, 10], the difference is that we suggest to calculate 
physical observables that are free from renormalza- 
tion uncertainties. Also, we demonstrate that the 
analysis of such correlation functions in continuum 
theory is aided by the conformal operator product 
expansion. 

The presentation is organized in follows. In Sec. II 



we introduce the coordinate-space correlation func- 
tion of two electromagnetic currents sandwiched be- 
tween vacuum and the pion. We calculate this cor- 
relation function to leading order (LO) in QCD per- 
turbation theory and discuss the physical interpre- 
tation. In Sec. Ill the state-of-the-art (next-to-next- 
to-leading order) calculation of this correlation func- 
tion is presented, including two-loop radiative cor- 
rections and nonperturbative twist-4 effects. The 
main result of this analysis is that the QCD cor- 
rections remain well under control for all pion mo- 
menta. In Sec. IV we discuss possible strategics and 
the potential accuracy of the extraction of the pion 
distribution amplitude from the coordinate space de- 
pendence of this correlation function, assuming that 
it is measured in a limited range of distances accessi- 
ble in lattice calculations. The corresponding setup 
and several possible generalizations are discussed in 
Sec. V, which also contains a summary and our final 
conclusions. 



II. EXAMPLE 

As an example, we consider the correlation func- 
tion of two electromagnetic currents sandwiched be- 
tween vacuum and the pion state 

T^, = {Q\T{j^{x)U-x)}n\p)) , (2.1) 

where 

2 1 - 

= -u{x)'yf,u{x) - -d{x)-f^d{x) , (2.2) 

which is the coordinate space analog of the pion 
transition form factor involving two photons. Note 
that one and the same correlation function (2.1) en- 
ters the pion decay tt^ — > 77 and the form factors 
7* — > 7r''7* for two virtual or 7* tt^j for one vir- 
tual and one real photons. Diff'erences arise at the 
stage when one goes over to the states with given 
momenta: depending on the virtuality, different co- 
ordinate space regions are emphasized/suppressed 
by the Fourier transform. 

At small space-like separations \x'^\ ^ 1/^qcd 
the amplitude in (2.1) can be calculated using the 
operator product expansion (OPE). To LO in the 
strong coupling, the answer is obtained from the di- 
agram in Fig. 1 and reads: 

= — -^/7re/ji/pcr^;^j3^T'(p • a;, 2; ) (2-3) 



with 




FIG. 1: Leading-order contribution to the correlation 
function (2.1) 



The relevant nonperturbative input is provided by 
the pion DA defined as 

{0\q{ain)jf,j5q{a2n)TT°{p)) = 

= iUp^ C due-^f"(«'^i+"'^= , (2.5) 



where /tt ~ 93 MeV is the pion decay constant, u = 
1 — u, Tip is a light-light vector, = 0, and /i is the 
normalization scale, which is set to be of the order 

In the following we write the pion DA as an ex- 
pansion over Gcgenbaucr polynomials, 

C5C 

Mu, fi) = 6uu Mf^)C^^^ (2u - 1) , (2.6) 

ri=0 

and use the fact that the Fourier transform 
8 / ciwwMe*''(2u-i)^3/2(2y _ 1) ^ 







= V^(n + l)(n + 2)zV-3/2j„+3/2(p) (2.7) 

leads to Bessel functions. To our accuracy we obtain 
then for the correlation function (2.4) 

Tip-x,x^)^-J2Ml^)^nip-x) (2.8) 



n=0 



with 



.f:„(p) =^"V27r 



T{p-x,x'^)= dwe'(2"-i)P-^(/.^(u,^). (2.4) 
Jo 



(2.9) 

We will view this result (2.8) as a partial wave ex- 
pansion of the correlation function T(p = p ■ x,x'^). 
(a group theoretical explanation will be given be- 
low). The partial waves J-n{p), expressed by Bessel 
functions Jn+3/2(p) with half integer index, are sim- 
ply given in terms of trigonometric functions, e.g.. 
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we have for n = 0: ^o(p) = 2[sin(p) — pcos{p))/p^. 
In particular this term, appearing as the leading 
one in the expansion (2.8), corresponds to the con- 
tribution of the asymptotic distribution amplitude, 
while T^ip) and J^i{p) take into account contribu- 
tions of the second, C^^^^'iu — 1), and the fourth, 
C4^^(2w — 1), Gegenbauer moments, respectively. 
Note that only even values of n contribute to the 
series because of C-parity conservation. The used 
approximation allows for a partonic interpretation 
of the correlation function T{p-x,x'^), which is obvi- 
ous: it corresponds to a probability amplitude of the 
valence quark distribution in the pion in the longitu- 
dinal distance ( "loffe-time" ) p = p ■ x at transverse 
distance x^. 

The partial waves Tq{p), — J-2(p) and T^ip) are 
plotted in Fig. 2. These oscillating functions are 




FIG. 2: From left to right: fo{p) (solid), -J^2{p) (dash- 
dotted), J-i(p) (dashed). 

strongly peaked at a position that moves with in- 
creasing index n to the r.h.s. and so the contribution 
of the asymptotic DA and the corrections arc well 
separated. On the other hand, large values oi p ■ x 
and, hence, of the pion momentum arc needed to 
probe higher-order Gegenbauer moments efficiently. 

To see this, we plot in Fig. 3 the correlation func- 
tion T(p, x^), given in the approximation (2.8), as 
a function of p = p ■ x for the asymptotic pion DA 
(dashed curve), and for the model with 02 = 0.25 
and two choices ip^ = 0.1 and (j)^ = —0.1, respec- 
tively, cf. Fig. 4. 

Note that the value of T{p, x^) at p = is equal to 
1/2 in our normalization (up to radiative and higher 
twist corrections discussed in the next Section) and 
does not depend on the shape of the pion DA. The 
deviation from the asymptotic DA (dashed curve) 
for a realistic model is significant starting p ^ 1. 
The dependence on 04, taken as 04 — 0.1 (dash- 
dotted curve) and 04 — —O.I (solid curve), becomes 
visible at distances p > 3 — 4. We remark that a 
more closer look, given in Sect. IV, reveals that a 
partial wave analysis might allow to access the forth 
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FIG. 3: The correlation function T{p,x^), leading-order 
(2.8), calculated using asymptotic pion DA (dashed), the 
model with 02 = 0.25 and two choices 04 = 0.1 (dash- 
dotted) and 04 = —0.1 (solid). 
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FIG. 4: Models for the pion DA as specified in Fig. 3. 

Gegenbauer moment even for p 3. With increasing 
p the correlation function changes sign and becomes 
negative, unless the fourth Gegenbauer moment is 
larger as the second one, which seems to be unlikely 
(not shown). 

The p-behavior of the models, plotted in Fig. 3, 
results from the cos-Fourier transform (2.4) of the 
DA (because for a symmetric DA the sin-component 
dies out) and can qualitatively be understood as fol- 
lows. The position of the first zero aX p = po gives 
us the effective width Tr/po of the DA in momentum 
space. For instance, a flat DA (p-niu) = 1 corre- 
sponds to po = TT, an equal momentum sharing DA 
07r('") ~ S{u — 1/2) leads to a constant correlation 
function, i.e., po oo, while the effective width of 
the asymptotic DA follows from pa w 1.437r. Com- 
pared to the asymptotic DA, for the realistic mod- 
els, shown by the solid and the dash-dotted curves, 
the first zero is shifted to the left, which indicates 
that these DAs are broader than the asymptotic one, 
cf. Fig. 4. The difference between the solid and the 
dash-dotted models, on the other hand, is mostly 
pronounced near the first minimum. The magni- 
tude of the correlation function at this minimum 
p = pi is a measure for the overlap of the DA with 
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the harmonic cos((2u— l)pi) which is maximized at 
ui = 1/2 — tt/pi ^ 0.2. It is large for the models 
with a two-hump DA like the one of Ref. [1] , where it 
was also argued that all higher Gegenbauer moments 
are negligible. It is much smaller for a dash-dotted 
model corresponding to a flat (albeit still oscillat- 
ing) DA as shown in Fig. 4. These oscillations are 
not significant and can be understood as an arti- 
fact of truncating the conformal spin partial wave 
expansion. In fact, a conjectured ADS/QCD model 
0ADS/QCD ^ (8/7r)^u(l -u) [11, 12] gives rise to a 
similar shape. We add to our discussion that very 
narrow DA's and those with nodes, leading to huge 
resonance effects, are clearly distinguishable in the 
correlation function from those shown in Fig. 3, even 
at smaller values of p. 



III. STATE-OF-THE-ART 

The previous discussion was mainly qualitative. 
In this section we present the state-of-the-art calcu- 
lation of the correlation hmction T(p- x, x^), defined 
in Eq. (2.3), in next-to-ncxt-to-leading (NNLO) of 
perturbative QCD and including the twist-four cor- 
rections. 

An adequate technical framework is provided by 
the conformal operator product expansion (OPE). 
Recall that the expansion in Gegenbauer polyno- 
mials in (2.6) alias the partial wave decomposition 
in (2.8) to LO are governed by conformal symme- 
try: the expansion is organized as the irreducible 
decomposition of the product of two currents into 
conformal operators. The partial waves specified in 
Eq. (2.9) are nothing but the LO Wilson-coefficients 
of the leading twist operators in the conformal OPE. 
They can be viewed as the Clebsch-Gordon coeffi- 
cients of the coUinear conformal group SL{2, R) and 
are labelled by the conformal spin j — n + 2, for 
a review see Ref. [14] . The advantage of such a de- 
composition is that the Gegenbauer moments do not 
mix under evolution in LO or, in other words, the 
conformal spin is a good quantum number to this 
accuracy. 

Taking into account radiative corrections to the 
leading twist Wilson-coefficients translates to the 
modification of the partial waves which become 
scale- (and scheme-) dependent: 

Tn{p) J^n{p, ~p'^x^;as{n)) ■ 

In turn, the higher-twist corrections yield x^ sup- 
pressed additive terms and for our present purposes 
it is convenient to reexpand them in terms of the 
leading-twist partial waves taken to leading-order. 



We end up with the expansion of the form 

T{p,x^)= (3.10) 
3 °° 

+ 1^' E ^n\i^^) Mp) + 0{x'a,) + Oix") , 

n=0 

where (jj^n^ are related to the matrix elements of 
twist- four operators, weighted with specific Wilson 
coefficients. Viewed this way, both the corrections 
to the partial waves and the twist-4 coefficients 
4>n\p'^) are specific for the considered correlation 
function. The Gegenbauer moments 4)n{p^) arc uni- 
versal, however, they depend on the scheme conven- 
tions. 

In what follows we consider separately perturba- 
tive and twist-four corrections in some detail. 



A. Perturbative corrections 

The ncxt-to-leading (NLO) perturbative correc- 
tions to the pion transition form factor were cal- 
culated in the minimal subtraction (MS) scheme in 
Refs. [15-17] and were supplemented by the evalu- 
ation of the logarithmical scale change [18-21]. In 
addition NNLO diagrams that are proportional to 
Uf, the number of light quark flavors, have been 
evaluated in the same scheme [22]. This result 
can be used to obtain the NNLO corrections that 
are proportional to /3o = (11/3)Ca — (2/3)n/, the 
first coefficient appearing in the QCD beta-function 
/3(ff)/.9 = -(as/47r2)/3o -I- 0{al). Finally, the con- 
straints imposed by the conformal symmetry, tested 
at NLO level [23, 24], have been used to obtain the 
missing terms at NNLO. A detailed NNLO analy- 
sis can be found in Ref. [26], see also Appendix A 
for a discussion of different factorization schemes. A 
Fourier transform of these results provides one with 
the radiative corrections to the correlation function 
(2.3) in position space. 

It is instructive to consider the results in the hy- 
pothetical conformal limit, in which the /3 function 
is vanishing (which means, technically, that the /3o 
proportional terms are omitted). In this case the 
modification of the partial waves, cf. Eq. (3.10), is 
entirely governed by conformal symmetry: 

T^^C^ia,) (_^V)*z"V^^^^±l^^(3.11) 

r(n + 5/2 + 7j2) /px-^ 
^ r(n + 5/2) I2J 
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As compared to the LO expressions (2.9), the p de- 
pendence is modified by the anomalous dimensions 
7„(as), which alter, e.g., the index of Bessel func- 
tions. In addition the normalization is changed by 
the factor 



Cnids) 



r(2 - 7„/2)r(i + 7- 

r(l + n + 7„/2) 



-c„(a,), (3.12) 



where c„(as) are the Wilson coefficients that ap- 
pear in polarized deep inelastic scattering structure 
function gi and are known to NNLO [27]. Since of 
conformal symmetry, a scale change does not lead 
to a mixing of conformal partial waves (3.12) or the 
Gegenbauer moments. 

One possibility to go beyond the conformal limit is 
to restore both the scale dependence of the coupling 
and the renormalization logs within the normaliza- 
tion condition of Ref. [26] . This scheme is discussed 
in Appendix A; we refer to it as the conformal sub- 
traction scheme, the CS scheme. We will use the 
renormalization group improved partial waves, de- 
fined in Eq. (A. 11), for the numerical analysis. In 
general, the anomalous dimensions govern the evo- 
lution of the Gegenbauer moments with respect to a 
scale change: 



-ln{as{^J))(j)n{^i^) 



(3.13) 



Even in the CS scheme which is designed to make 
maximal use of the conformal symmetry, we ex- 
pect that the Gegenbauer moments will mix un- 
der evolution, indicated by the ellipsis. This mix- 
ing is induced by the trace anomaly and gives rise 
to a (as/27r)^/3o ln(— a;^/x§) proportional contribu- 
tion which, however, vanishes at a reference scale 
—x^ = l/no- We expect that this unknown NNLO 
mixing effect is small and can safely be ignored. 

The partial waves in the conventional MS scheme 
can be obtained from those in CS by the appropri- 
ate transformation. In particular the Gegenbauer 
moments in these two schemes are related as 



/CS 



/MS 



n-2 

as \ ^ 

27r ^ 

Tn=0 



Sill0™' + O(a^). (3.14) 



The matrix i?im is known explicitly and is given in 
Eq. (A. 16) of Appendix A. This relation is also valid 
in the hypothetical conformal limit. Note that the 
MS Gegenbauer moments are given by a finite sum of 
those in the CS. Consequently, plugging Eq. (3.14) 
into Eq. (3.10), one sees that in the MS scheme all 
conformal partial waves (3.11) are excited even if a 
truncated model for the DA is used. We stress that 
the scheme independence of the correlation function 



is only guaranteed if the Gegenbauer moments are 
rotated at the given input scale. Taking the same 
model in different schemes will in general lead to 
different predictions for the correlation function. 

As we will see in Sect. IV, such scheme-dependent 
mixing effects are numerically important. For illus- 
tration, consider the asymptotic DA which we define 
here as the zero mode of the evolution kernel in a 
given scheme and order of perturbation theory. By 
construction, the asymptotic DA does not evolve un- 
der scale transformations. In the CS scheme to NLO 
this amounts to the same choice of the Gegenbauer 
coefficients as in LO: ipo = 1 and = for n > 0. 
The one-loop corrections to the correlation function 
(3.10) are then given by 



1 - 



0{al) 



Mr)- (3.15) 



The current conservation implies 70 (a^) — 0, pro- 
tecting the p-dependence of the lowest partial wave 
from radiative corrections. In the MS scheme, on the 
other hand, the Gegenbauer moments of the asymp- 
totic DA for n > 2 arc nonzero and arc obtained by 
the rotation (3.14): 



-as-MS, 



a^(/i) 8Ci.(2n-h3) 
27r n(ri-t- l)(n-|-2)(n + 3) 



(3.16) 



X Si{n + 2) 



2{n + 1) 



whereas the normalization, i.e., (pQ = 1, re- 
mains unchanged. The excitation of higher Gegen- 
bauer moments is, however, rapidly decreasing 
with n: For n = {2,4,6,8} we find the values 
{0.124,0.039,0.018,0.009} in units of as{p), respec- 
tively. A more extensive discussion about scheme 
dependence can be found in Ref. [29]. 

Finally, the limit p ^ of the correlation function 
(2.3) yields the analog of the Bjorken sum rule in 
polarized deep inelastic scattering: 



1 



(3.17) 



In this limit, which is independent on the factoriza- 
tion scheme, any information contained in the DA is 
washed out. Hence, this sum rule is a pure QCD pre- 
diction, which is known as pcrturbative expansion to 
three loop order [28]. 

After these general remarks, we consider the per- 
turbative corrections for the particular models of the 
pion DA specified in Sect. II. We assume that the 
models are defined at the scale 1 GeV and set the fac- 
torization scale to p? = — l/x^. For the demonstra- 
tion we consider two current separations: a larger 
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one -x^ = 1/GeV^ [«(0.2fm)2] and a smaller 
one -x^ = l/4GcV^ [w (0.1 fm)^] . In the lat- 
ter case evolution of Gcgcnbauer moments is taken 
into account: The leading logs arc rcsummcd and 
non-leading ones are consistently combined with the 
as /27r power expansion of the correlation function 
(3.10), for details see, e.g., Ref. [26]. The couphng 
is specified at the normalization point ^ = 2 GeV to 
be as = 0.36 and as = 0.3 at LO and NLO, respec- 
tively. 
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FIG. 5: NLO predictions for the correlation function 
(2.3) in the MS scheme (thin) and rotated to the CS 
one (thick) at tlie scale —x^ — 1/GeV^ (lower panel) 
—x^ = 1/4 GeV^ (upper panel) compared to the LO re- 
sults (dots) for the three models of the pion DA shown 
in Fig. 4. 

The NLO results are shown in Fig. 5 for —x^ = 
1/GeV^ and -x"^ = l/4GeV^ in the lower and up- 
per panel, respectively. Notice that the pertur- 
bative corrections lead to a universal shift of the 
value at /5 = 0, which is 1/2 in LO, downwards to 
1/2(1 — as/ir) which is about -15% for the larger 
and -9% for the smaller separation between the cur- 
rents, respectively. This effect is entirely due to the 
negative NLO correction in the sum rule (3.17). 

Besides the overall normalization, perturbative 
corrections result in a slight change of the p depen- 
dence, seen as a shift of zeros and extrema to larger p 
values (compare thin and dotted lines) . This scheme 
dependent effect arises from the excitation of higher 
conformal partial waves, as discussed above, and is 
more pronounced at larger separations (lower panel). 
Note that the p-dependence corresponding to the 



asymptotic DA (dashed curves), evaluated in the CS 
scheme (thick curves), is not affected, see Eq. (3.15). 
For the other two models in the CS scheme only a 
slight shift of extrema appears, which is visible in the 
lower panel at larger values of p (compare thick and 
dotted lines). A more quantitative look reveals that 
in general the absolute size of radiative correction 
grows somewhat with increasing p. 

Note that at smaller distances the model depen- 
dence slightly weakens, compare the magnitude of 
solid or dash-dotted curves at the first minima in the 
upper and lower panel. This is caused by evolution, 
which does not change the value at p = 0, however, 
reduces Gegenbauer moments with a strength that 
grows with increasing conformal spin. Obviously, 
even at —x'^ = l/4GeV^ the model predictions re- 
main clearly distinguishable. 

The NNLO corrections evaluated in the CS 
scheme (not shown) lead to a further decrease of 
the normalization, so that the net reduction for 
p = is, compared to LO, —25% and —13% for 
-x"^ = 1/GeV^ and -x^ = l/4GeV^, respectively. 
The modification of the p dependence is, however, 
negligible as compared to the NLO result in the CS 
scheme, and is hardly visible even for the larger sepa- 
ration —x^ = 1/GeV^. As mentioned above, mixing 
effects due to the NNLO evolution are not taken into 
account, but are expected to be tiny. 

B. Higher- twist Corrections 

Higher-twist corrections are in general suppressed 
by full powers of the separation x'^ between the cur- 
rents. In particular, including twist -four contribu- 
tions one obtains 

T{p-x,x^)^ [ due*(2"-i)P-^t(M,.T2), (3.18) 
Jo 

where 

t{u,x'^)=(j)Tr{u) + — 04(m) — .T^ / dai / da2 
4 Jo Jo 

X — <^4{a) + - y -<^4{a) . 

aa a% 

(3.19) 

The leading-twist perturbative corrections arc not 
shown for brevity. The notations correspond to 
Ref. [30], where one can find the definitions and 
the corresponding expressions for the twist-four dis- 
tribution amplitudes to the next-to-leading confor- 
mal spin accuracy. In the last term ai,ck2 and 
as = 1 — ai — q;2 are the quark, antiquark and gluon 
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momentum fractions, respectively. Including the 
contributions of the lowest and the next-to-lowest 
conformal spin one obtains [30, 31] 

(j)4{u)~—^5'^u^u^ + 21(5^tJ47r|uu(2 + 13mm) 

+ [2u^{6u'^ - 15m -f 10) In m] + [u ^ m]| , 



I 21 , 
21 



*4(a)=120Q;iQ;2Q;3(Q;i - q;2)i5^— t^47r 



(3.20) 



The nonperturbative parameters 5"^ and 0^4^ are de- 
fined as reduced matrix elements of local operators, 
for example 



(3.21) 



Numerical estimates for these matrix elements are 
available from QCD sum rules: 



Si = (0.18±0.06) GeV^ 



LUi^ = 0.2±0.1 (3.22) 



at the scale 1 GeV [S^ = (0.14 ± 0.05) GeV^ uji^ = 
0.13 ±0.07 at the scale 2 GeV], see Ref. [30] and the 
references therein. 

Note that beyond the leading conformal spin ac- 
curacy the twist-four contributions are not polyno- 
mials [31]. Plugging Eq. (3.20) into Eq. (3.19), we 
find, however, that the next-to-leading conformal 
spin contributions completely cancel each other and 
the leading conformal spin yields the total O(x^) 
correction: 
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2„-:2 



t(M,.T^) = (j)2iu)+X^-di ■ 30m^M 



(3.23) 



This can be effectively rewritten as the correction to 
the first two leading twist Gegenbauer coefficients 
1^1 + (.S/9)x'6l 02 ^ </,2 - (8/9)(l/6)x2<52. 
Hence, we have 

c|^i'\^^) = lsl{^,), 4^\^,)^^-^sl{^^). (3.24) 

The upper bound for the contributions to (3.20) of 
higher conformal spins can be obtained using the 
renormalon model of Ref. [37] . These extra contri- 
butions mainly influence the large p-behavior of the 
correlation function while the numbers in (3.24) are 
not affected. We see that, to LO accuracy, higher- 
twist correction produce an additive shift in the 
physical values of the Gegenbauer moments of the 
pion DA, which is calculable, at least in principle. 
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FIG. 6: LO predictions for the correlation function (2.3) 
with (thick) and without (thin curves) twist-four correc- 
tions for two choices of the quark-interquark separation; 
-x^ = 1/GeV^ (lower panel) and -x'^ = 1/4 GeV^ (up- 
per panel). Models for pion DA are the same as used for 
the LO predictions, shown by dotted curves, in Fig. 3. 



Beyond LO, one may define effective partial 
waves. 



(l>n{p) 



(3.25) 

that additivcly combine pcrturbativc and higher 
twist corrections together. Here, the evolution of the 
twist-four coefficient is governed by the difference in 
the LO anomalous dimensions of the corresponding 
operators: 



<l>n{p) 



4>niPo) ' 



OisiPo] 



(3.26) 



where 7q°^ = and 72"'' = 50/9. Note that these ef- 
fective partial waves depend on the non-perturbative 
quantities, which induces some model dependence. 

As it is demonstrated in Fig. 6, twist-four con- 
tributions are significant at —x^ = 1/GeV^ (lower 



panel) and much less important at 



l/4GcV^ 



(upper panel). In particular, for —x^ = 1/GcV^ 
the value at p = decreases by about 16% which 
is comparable to the NLO perturbative correction, 
whereas for — ~ l/4GeV^ the decrease is only by 
~ 3.5% which is roughly one third of the respective 
NLO effect. Note that the p dependence of separate 
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partial waves is not affected by the twist-four con- 
tributions. Their magnitude is changing, however: 
the second partial wave is enhanced and the lowest 
one somewhat suppressed by the twist-four correc- 
tions. These effects depend linearly on and are 
additionally suppressed by evolution. 



> 
O 



> 
O 




FIG. 7: Combined NNLO and twist-four predictions for 
tire correlation function (2.3) in the CS scheme (thick 
curves) compared to the LO results (thin curves). The 
factorization scale is set to the quark-interquark sepa- 
ration, ii^ = — which is —x^ = 1/GeV^ {—x^ = 
1/4 GeV^) for the lower (upper) panel, respectively. 
Models for pion DA are the same as in Fig. 3. 

Our final predictions combining the perturbative 
effects to NNLO accuracy in the CS scheme and the 
twist-four corrections are shown by thick curves in 
Fig. 7 and compared to the corresponding leading- 
order leading twist predictions (thin curves). For the 
separation ~ l/4GeV^ the overall correction is 
rather small. The normalization at p = is reduced 
by about 17%, whereas the zeros and positions of 
extrema are essentially fixed. We remind that there 
will be small flow of these points in the MS scheme, 
cf. upper panel in Fig. 5. The overall corrections for 
—x^ = 1/GeV^ arc much larger. 

The conclusion is that the correlation function 
(2.3) can be calculated in QCD to a high accuracy. 
In particular for — < l/4GeV^ ~ (0.1 fm^ the per- 
turbation theory works rather well and the higher- 
twist effects do not exceed 3%. 

The next question to address is whether the 
knowledge of the correlation function (2.3) can be 
used to constrain the shape of the pion DA and, for 



example, allow for an accurate determination of the 
few lowest Gegenbauer moments. 



IV. REVEALING THE PION DA 

Let us first consider the QCD prediction for the 
normalization of the correlation function in the 
vicinity of p = 0. Here only the lowest partial wave 
contributes, while higher ones are suppressed as 



\J'n{p)\ < 



(n+ l)(n + 2)0F 
4r(n + 5/2) 



(4.27) 



Note the factorial suppression factor. For the fa- 
vored value of (j)2 ~ 0.25, we find for p ~ 1/4 [~ 1/2] 
that the contamination of the second partial wave is 
already below three per mil [about 1%], respectively. 
All higher partial waves can certainly be omitted. 
We remind that a strong suppression of higher par- 
tial waves occurs also in momentum space, if both 
photons have a large virtuality [25, 26]. It follows 
that to three per mil [1%] accuracy the sum rule 
(3.17) can be extended for p < l/4[l/2]: 



+0iasx^) + Oix^) , 



(4.28) 



where {3/A)To{p) = l/2--pV20+O(/). We empha- 
size again that in the CS scheme the radiative cor- 
rections to the lowest partial wave are absent at the 
normalization point. We expect that the excitation 
of higher partial waves due to the evolution yields a 
neghgible 0{p'^a'^) effect, see also Eq. (3.14) and the 
discussion of the MS mixing. Hence, in this scheme, 
the perturbative corrections can be borrowed from 
the Bjorkcn sum rule, known to order [28]. For 
shortness wc display the two-loop result: 

Cf(V-^a.(M)) = 1-^-^^(4.29) 



Note that due to the Fourier transform the argu- 
ment of the rcnormalization logs is decorated with a 
transcendent number 6^^"=^^^ « 1.166994- • •. 

To predict the normalization of the correlation 
function, we have to specify the QCD coupling. The 
world average value of ag , given at the Z-boson mass 
scale, reads to three loop accuracy [32] 

ai^\fi = 91.18 GcV) = 0.1184 ±0.003. (4.30) 
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FIG. 8: Normalization uncertainty of the correlation 
function (2.3) for p — 1/4. The upper error band entirely 
arises by the uncertainty of the running coupling, while 
the lower one takes also into account twist-four contam- 
ination within the central value, given in Eq. (3.22). 



In the backward evolution to a lower scale we take 
into account quark thresholds, treated in the stan- 
dard way [33]. For p ~ and the current separation 
— ~ l/4GeV^ ^ (0.1 fm)-^ we obtain, including 
three-loop and higher-twist effects:: 

T(0, -l/4GeV2) = 0.4301° ™? - 0.0181° °°^ 

(4.31) 

where the number of quarks is set to n/ = 4. The 
predicted twist-four correction is comparable to the 
error induced by the uncertainty of a^, sec Fig. 8, 
where the dependence of T(0,a::^), including the er- 
ror bands, is displayed as a function of x^. 

Next we discuss the access to the Gegenbauer mo- 
ments. We assume that the correlation function is 
measured in a certain interval of pion momenta alias 
an interval p, po < p < Pmax- Since an overall nor- 
malization might be difficult, we consider the ratio 



R{po,p,x^) 



T{po,x'^) ' 



Po<p, (4.32) 



which is normalized to one at p = pQ. Requiring that 
po < 1 allows to truncate the partial wave expansion 
in the denominator at the second term ^ 4'2- As the 
result, R{pq, p,x^) is essentially a linear function of 
(j)n with n > 2. 

The maximal accessible value of p, given by Pmaxj 
is limited by the pion momentum that is feasible in 
a lattice calculation. We assume here that pmax < 3, 
i.e., lower than the first zero of the lowest conformal 
partial wave, cf. Fig. 2. Staying away from the zero 
allows one to consider the subtracted ratio to access 
directly the second Gegenbauer moment: 



Mipo,p,x^) = 1 - 



R{po,p,x^] 
T^o{pa,p) 



Po < P < 3 , 

(4.33) 



where 7^„(po,p) = ^nip) / ^o{po) ■ Since 7^o(po,p) 
only decreases by at most a factor three in the whole 
range p < 3, this weight does not introduce any 
numerical instability. The conformal spin expansion 
of A4 in terms of ratios of the effective partial waves 
(3.25) reads: 



M 



^ [7^f (po)-7^f(p)] 02 _ 7^f (p)(/>4 



l+Tlfipo)(b2 



where 



l + (PO)02 ' 

(4.34) 



nfip,xy- 



•^r(p,^^) 
:Ffip,x^) 



(4.35) 



Here we neglected all contributions in (f>Q, 08, etc. 
from the expansion in the numerator. The corre- 
sponding truncation error is estimated to be on the 
per mil [percentage] level for p < 2 [3], respectively. 
For example, for p = 3 we find to LO accuracy: 

M - M'-'^^ - 3 • 10^-^04 + 7 • 10-^06 + 3 • 10-^08 • 

(4.36) 

Here the first term on the r.h.s. comes from the trun- 
cation in the dominator (assuming po = 1/2) and 
can safely be neglected, same as the contribution of 
08- The remaining term in (f>e effectively induces a 
small additive uncertainty in determination of the 
first two Gegenbauer moments. For instance, as- 
suming values of A4 are known for two different 
values p = 2 and p = 3, we find by solving the set of 
two linear equations: 

S(j)2 ~ -0.0106 S(j)4 w -0.1306 ■ (4.37) 

Assuming 06 ~ 0.1 the errors induced by the trun- 
cation of the partial wave expansion are in this ver- 
sion negligible for 02 and probably on the one per- 
cent level for 04. We remind that convergency of 
the Gegenbauer expansion in the vicinity of u = 
{u — 1), which is required if the DA vanishes at the 
end points, leads to an upper bound for the large- 
n behavior of the coefficients |0„| < const /n^ with 
p > 1, so that large values of 0„ with n > 4 are 
increasingly unlikely. Clearly, a fitting procedure in- 
cluding 06 is also possible. This is demonstrated in 
Fig. 9, where the effect on Ai is shown for the 06 
varied in the interval —0.2 • • • 0.2, which we consider 
as an overestimation. 

Radiative corrections to the ratio (4.32) of the cor- 
relation function are much milder as compared to 
their overall normalization. A precision analysis re- 
quires a careful specification of the running coupling. 
We take the world average value (4.30) and calculate 
the corresponding value at the scale p, = 2 GeV by 
the tree loop evolution equation within Uf =4: 



as{p = 2GeV)=0.304 



+0.024 
0.021- 



(4.38) 
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p = p ■ X 

FIG. 9: Subtracted ratio (4.33) of the correlation func- 
tion (2.3) calculated at LO as a function of p with 
po = 1/2 and -x^ = 1/GeV^. Models of the pion DA 
are same as in Fig. 3, where the splitting of curves shows 
the possible (overestimated) contribution of higher con- 
formal partial waves with n > 4. 




0.5 1 1.5 2 2.5 3 



p = p ■ X 

FIG . 10; Same as in Fig. 9 to NLO for the MS (thin) and 
CS radiative (thick) schemes. The variation within the 
predicted twist-four corrections are shown as error band. 
Here —x^ = 1/4 GeV'^ and the coupling is specified by 
the central values (4.38). 



We note that the resulting central value and errors 
of Us at two loop accuracy are compatible with the 
given one; it is, however, in the heart of perturba- 
tion theory that the specification of Us to LO ac- 
curacy goes hand in hand with large uncertainties. 
As discussed in Sect. Ill A, the asymptotic DA at 
NLO in the MS differs from that in the CS scheme, 
mainly for the n = 2 and n = 4 moments, see 
Eq. (3.16). This effect is seen in Fig. 10, where 
the subtracted ratio Ai (4.34) calculated in the MS 
scheme using the LO asymptotic DA (thin dashed 
curve) clearly deviates from zero. This deviation 
is removed for the corrected asymptotic DA (3.16), 
for which the n = 2 moment is « 0.04. For the 
same reason, the other two model predictions in the 
MS scheme (thin curves) are systematically shifted 
downwards compared to the CS (thick curves). A 



small deviation from zero for the prediction corre- 
sponding to the asymptotic DA in the CS scheme 
(thick dashed curve) comes from the excitation of 
the n = 2 partial wave by a twist-four contribu- 
tion. The shown narrow error bands arise from the 
variation of the twist-four parameter in the range 
5l{^l = 2GeV2) = 0.09. .. 0.19 GeV^. 

Since perturbative corrections to the partial waves 
in the CS scheme only induce a small modification 
of the p dependence in the considered range (absent 
for the lowest one and on the one percent or lower 
level for higher ones), the extraction of the first few 
Gcgenbaucr moments can be considerably simpli- 
fied. To this end, we define the effective Gcgenbaucr 
moments: 

^off ^ C'„(((p))„|a,,-a;V^)0„ 

Co(a.,-xV') + (8/9)^252 ' ^ 

which depend on the average value of p in the given 
range {{{p))n ^ 2 for p^ax ^ 3). 

The subtracted ratio (4.34) in the range po < p < 
Pmax can be written to a good accuracy in terms of 
the effective moments as 

_^(4) ^ [7^2(po)-7^2(p)] j>f _ 7^4(p)<l>f 

l + 7^2(po)$f l + 7^2(po)$f' 

(4.40) 

where TZ„ {p) = Tn{p) / J^o[p) is the ratio of the LO 
partial waves (2.9). 

The relation between the effective moments and 
the physical Gegenbauer moments in the expansion 
of the pion DA includes radiative and higher twist 
corrections. As an illustration, for ((p))„ = 2.2 and 
ai p = 2 GeV we find 

02(2GeV) = [O.96ltE5:;5S]^ _ o.03t°:°}] $f 

+ 0.005t°:l]°2 + . . . , (4.41) 

04(2 GeV) = [0.902l°;r - 0.031^1] < + •••• 

The first term in the square brackets corresponds 
to the perturbative correction while the second one 
arises from the twist-four contribution. Note that 
the normalization change due to radiative correc- 
tions is only on the level of 4% for (t)2 and 10% for (/)4, 
respectively, and is comparable with the higher twist 
correction. In NNLO we find a slight increase of the 
normalization 0.961 0.982 and 0.902 0.951 (for 
Uf = 4), whereas the error is drastically reduced. 
The residual factorization and the renormalization 
scale dependence also prove to be very small. 

As the final step, the Gegenbauer moments in the 
MS scheme to NLO follow from 

0r=02 + ^^+O(a^), (4.42) 
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V. CONCLUSION 

A typical setup for a lattice calculation of the cor- 
relation function (2.1) would be to take the sepa- 
ration between the two currents purely spacelike, 
—f {0, a?}, and integrate over the cm. position 
of the currents with a certain three-dimensional mo- 
mentum. Since the pion coupling to the source is 
usually momentum dependent, the simplest strategy 
would probably be to consider ratios of the correla- 
tion function calculated for the same pion momen- 
tum but different x (and different cm. momenta of 
the pair of currents), thus creating the set of "data 
points" with correlated values of x^ = — |xp and 
p = p ■ X. In this way the dependence on the pion 
coupling cancels out. The crucial condition for such 
a calculation is to have a sufficient "lever arm" in 
p. For a lattice spacing a the momentum cannot 
be larger than 

\p\ < , 

which translates to the restriction on the maximal 
accessible value of p 

Pmax ^ 

a 

In other worlds, p cannot be larger than (half of) 
the separation between the currents in lattice units. 
Since, on the other hand, pcrturbative treatment of 
the correlation function is possible for \x\ < 0.2 fm, 
this requirement translates to a necessity to work on 
a fine lattice with a < 0.03 — 0.05 fm. This presents 
a considerable challenge for an unqucnched calcula- 
tion. The advantage of this formulation is, however, 
that the problem of (nonperturbative) renormaliza- 
tion of lattice operators is avoided altogether. 

For a demonstration, we have chosen in this work 
to consider the correlation function of two electro- 
magnetic current in which case the pcrturbative ex- 
pansion is known to the NNLO accuracy. This 
choice may be inconvenient for a lattice calculation 
because of the epsilon-tensor appearing as a prefac- 
tor in (2.3). This problem can easily be avoided by 



considering a correlation function of a vector and an 
axial- vector current, e.g., 

- (0|r{g(x)7^g(a:)g(-x)7,75g(-x)}^0(p)) . 

(5.43) 

More importantly, the same strategy can be used 
to extract information about the B-meson distribu- 
tion amplitude in the heavy quark limit from the 
lattice measurement of the correlation function of 
the type 

{0\T{K{xh^q{x)q{-x)Tq{~x)}\B^,{p)) , (5.44) 

where v is the heavy quark velocity and hy an effec- 
tive heavy-quark field operator, which can be calcu- 
lated in terms of the B-meson DA using soft-collincar 
effective theory [38-40]. The standard method (cal- 
culation of the moments) is not applicable in this 
case, since the existing definition of the B-meson 
DA relies entirely on a pcrturbative factorization; 
the relation with the Wilson operator expansion is 
lost unless an additional energy cutoff is introduced, 
see [41, 42]. 

Finally, nucleon distribution amplitudes [43, 44] 
can be studied in the same manner, from the corre- 
lation functions involving a local baryon current. 
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APPENDIX A: PERTURBATIVE 
EXPANSION 



The radiative corrections in the position space might be obtained from those in the momentum space by 
Fourier transform: 



^-^T{p-x,x)=Jj^^e^ Tiu;,-q), (A.l) 
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where the hard scattering amphtudc T{uj, ~q^) is expressed by the asymmetry variable uj = —q ■ p/q'^ and 
q ~ {qi — q2)/'^- Note that to LO accuracy the normaUzation is chosen to be T{uj = 0) = 1. The correlation 
function is given by the partial wave expansion: 



(A.2) 



In the hypothetical conformal limit the running coupling is replaced by its non-trivial fixed point value, at 
which the /3-function vanishes. Then the predictive power of the conformal operator product expansion can 
be used to obtain the partial waves: 



Tn[uJ, ~^l /q |7„) = Cn{a,)——r— — , — 



2"+2r(7i + 5/2) y-q"^ 



(n + l + e)/2,(n + 2 + e)/2 
(n + 5 + e)/2 



Here 



7n(as) 



-In 



(0) 



27r ■ (27r)2 
are the anomalous dimensions and 

.(1) _ 



-3 



1 



27r " 



+ 



Cf Sfin + 1) - S2in + 1) + -Siin + 2) + 



(n+ l)(n + 2) 

3-251(71) 



'-=l,^/■2 

(A.3) 



(A.4) 



2(n+ l)(n + 2) 



- t; (A.5) 



are the Wilson coefficients, normalized to one (cn'' = 1) at LO, of the polarized deeply inelastic structure 
function 171 for the flavor non-singlet sector. Both of them are evaluated in the momentum space within 
the MS scheme. The NLO anomalous dimensions can be found, e.g., in Ref. [34], and the few lowest one to 
NNLO in Ref. [35]. The NNLO expressions (or numerical values for the few lowest) Wilson-coefficients c„ 
are obtained from Ref. [36]. All numbers, used in this paper, are collected in Ref. [26]. 

To transform the conformal predictions (A.3) into the position space, we first utilize the standard integral 
representation for hypergeomctric functions and in addition a quadratic transformation: 



(n+l + e)/2,(7i + 2 + e)/2| 
(n-t-5 + e)/2 ' 



r(2n- 



2e) 



r(n-t-2 + e)2 



du 



[(1 - u)u 



n+l+e 



(1 + UJ — 2uju) 



n+l+e 



(A.6) 



Plugging Eq. (A.3) into Eq. (A.2), and interchanging u and q integration in Eq. (A.l), the Fourier transform 
can now be easily performed by means of: 



(27r)2 



1 



r(2-e) Xp{-x^fi^Y 



167r2r(l + n + e) {-x^)^ 



(-ip-a;)"e'(2«-i)^-P, (A.7) 



where q ^ q + (2u— l)p/2. Here and pp terms are neglected to leading twist accuracy (the latter drops out 
after contraction with the Levi-Civita tensor). Finally, the u-integration leads to the partial waves (3.11) in 
position space within the Wilson coefficients (3.12). 

The conformal symmetry breaking due to the trace anomaly appears the first time as aj.Po correction. 
It contains the renormalization logs, needed to restore the renormalization group invariancc in the pcrtur- 
bativc expansion with respect to a running coupling as{p), and induces a mixing of Gcgenbaucr moments, 
which depends on the factorization scheme. It is purely conventional how we are dealing with this mixing. 
Appealing possibilities are 

i. combining the conformal predictions (A.3) with the diagrammatical MS result [22], 

ii. the partial waves matches at the normalization point —q^ = /i^ the conformal predictions (A.3), 



iii. the Gcgenbaucr moments evolve autonomously. 
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Certainly, the latter case is the most appealing ones, however, it requires the full NNLO result for the 
trace anomaly induced term (in an arbitrary scheme). Unfortunately, the corresponding piece of the NNLO 
anomalous dimensions is unknown. Hence, we are only able to follow the first two suggestions, where the 
mixing due to evolution remains unknown [26] . It turned out that the numerical differences of the constant 
(3Qa1l2 proportional terms of these both schemes, are moderate. Moreover, it is expected that the mixing 
effect due to the evolution is tiny. We remind that in the MS scheme such a mixing under evolution already 
appears at NLO and leads in fact only to small numerical effects. 

Following the second suggestion, we present here the NNLO corrections in the so-called CS scheme. To 
do so we restore in the conformal predictions (A. 3) the scale dependence of the coupling, i.e., as as^J-), 
and the renormalization logs, which are governed by the renormalization group equation 



9 ar ^ d 



dg 



TniuJ-.-n'^ /q^\as{^)) = 7„(as(/i))J^„(cJ, -y? / q^;as{^.)) + 0{al) , 



(A., 



approximated to NNLO. The remaining freedom is fixed by the requirement that at the normalization point 
—q^ = p?' the partial waves coincide with those in Eq. (A. 3). A more detailed discussion is given in Ref. 
[26]. The so found result. 



(A.9) 



is understood as a perturbative expansion up to NNLO accuracy, where the addenda, 



A(^'«(c.,-MVg2) = ln(-MV'?') 



d 



(A.io) 



as=0 



restores the renormalization group invariance. 

Again utilizing Eq. (A. 7) as expansion with respect to e, the Fourier transform of the addenda (A.IO) 
is straightforward. Combining the resulting expression with the perturbative expansion of the conformal 
predictions (3.11) and (3.12) we find up to order a^: 



Mp), (A.ll) 



where Tn{p) is the LO partial wave (2.9) and 



The 'shift' functions 



sMfo -x^a^) - ( 2 2^. r(2-.)r(l + n).F„(p,-xV|26) 



(A.13) 



c=0 



depend on the variables —x^p? and p. To obtain the renormalization group improved Wilson-coefficients we 
restore the scale dependence of the coupling, i.e., as as{p), and take into account 



C(2,/3) 



^Ln„(-.T2M') - C(i)(p,-a;V')) Ln„(-xV') + (1 - S^in)) 



(0) 



(A.14) 



where Ln„(— x^/j,^) = ln(— x^^^) — S'i(n) + 27E — 1. We finally remark that the following substitutions restores 
the result in momentum space: 



(" + !)(» + 2)V^ .„ ^({n + l)/2, {n + 2)/2 



2"+2r(?i + 5/2) 



(n + 5)/2 



UJ 
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(ri+2c+l)/2,(n+2£+2)/2 I 
n+2e+5/2 I 



i+2)/2 I 



(A.15) 



p / (n+l)/2,(n^ 
2^1 n+5/2 ry e=0 

The NLO corrections to the correlation function in the MS scheme might be obtained by Fomier transform 
from the known result in momentum space [15-17]. Representing the hard-scattering part as a convolution 
of the LO hard-scattering part within some kernels, see Ref. [23], the transformation can be straightforwardly 
performed by utilizing Eq. (A. 7) for n = 0. Alternatively, we can simply rotate the partial waves (A. 11) in 
the CS scheme, expanded up to NLO. Note that the trace anomaly does not enter in this approximation. 
The non-vanishing entries of the rotation matrix, used in Eq. (3.14), 



(1) - 



2(2n+3)Cf 
(n+l)(7i + 2) 



(m+l)(TO+2) 



lAr, 



1' 



(0) 



2Cp 



+ Anm - Sl{n+ 1) 



(n — m)(?i+m + 3) 
appear for n > m and n — m-even, where 

An,n ^Si{{n + m + 2)12) - Si[{n - m - 2)/2) + 2Si{n - m - 1) - Si{n + 1) 



(A.16) 



(A.17) 
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